Experimental determination of the absorption strength in absorbing chaotic cavities 



O 

o 

(N 
•*-> 

o 
O 



G. Baez, 1 M. Martmez-Mares, 2 and R. A. Mendez-Sanchez 3 

'Area de Fisica Tedrica y Materia Condensada, Universidad Autonoma 
Metropolitana-Azcapotzalco, A. P. 21-267, 04000 Mexico D. F., Mexico 
2 Departarnento de Fisica, Universidad Autonoma Metropolitana-Iztapalapa, A. P. 55-534, 09340 Mexico D. F., Mexico 
1 Instituto de Ciencias Fisicas, Universidad Nacional Autonoma de Mexico, A. P. 48-3, 62210 Cuernavaca Mor., Mexico 

(Dated: February 2, 2008) 

Due to the experimental necessity we present a formula to determine the absorption strength by 
power losses inside a chaotic system (cavities, graphs, acoustic resonators, etc) when the antenna 
coupling, always present in experimental measurements, is taken into account. This is done by 
calculating the average of the absorption coefficient as a function of the absorption strength and the 
coupling of the antenna to the system, in the one channel case. 

PACS numbers: 73.23.-b, 03.65.Nk, 42.25.Bs,47.52.+j 
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I. INTRODUCTION 

A great interest on chaotic systems with internal losses, 
or absorption, is due to the fact that they are good can- 
didates to verify Random Matrix Theory (RMT). The 
effect of losses has given rise to many theoretical and ex- 
perimental investigations on transportproperties (for a 
review see Ref. [lj, see also @,i SSI). In particular 
microwave B, Q and acoustic resonators 0] , microwave 
networks |lOj and elastic systems [ll[ are excellent tools 
for that purpuse. 

One of the most important things in absorbing systems 
is to quantify the degree of the losses suffered in an ex- 
perimental situation. Although in some experiments the 
absorption strength 7 can be p artially controlled by at- 
taching aditional antennas [13] or introducing absorber 
materials inside the system [2|, the intrinsic absorption 
is not directly under experimental control. However, it 
can be tunned to fit the experimental data [H, [13, [H[ . 
In a single-mode port microwave resonator setup of Refs. 
El 1 the one channel case, 7 was determined as the 
one that reproduces the experimental value of the average 
of the reflection coefficient (R), which is a monotonically 
decreasingfunction of 7. A similar procedure was used 



in Ref. b| but with the transmission coefficient T in 
a two single-mode ports. In that case, 7 was chosen as 
the value in which the theoretical distribution of T fits 
the experimental one. Also, the value of 7 was estimated 
from the fidelity [l6| . 

The imperfect coupling of the antennas to the cavity 
has a similar effect as the absorption and unfortunately 
the last procedure can not discern what part of 7 is due 
to losses and what to the rejected part of the wave that 
never enter the cavity. In this sense an imperfect coupling 
mimics absorption and a procedure that gives the correct 
value of 7 is needed. Here we present a semi-analytical 
formula to calculate 7 which takes into account the cou- 
pling for the one channel case. This is helpful in exper- 
iments with microwave networks and one port cavities. 
Although the procedure is the same for higher number 
of channels, it is not possible to give a result due to the 




FIG. f : Sketch of a flat chaotic cavity. So describes the scat- 
tering of the cavity with perfect coupling of the antenna repre- 
sented as a flat waveguide. S describes the scattering through 
the system cavity plus a barrier that model the imperfect cou- 
pling. 



absence of knowlegde of the probability distribution of 
the scattering matrix 1] . 

In the next section we summarize the existing theory 
used to describe the scattering process through chaotic 
cavities with losses, taking the imperfect coupling to a 
single- mode port into account. Sect. IIIII is devoted to 
calculate the absorption strength 7 in terms of the av- 
erage of R the coupling intensity. We conclude in Sect. 



II. THE S-MATRIX AND ITS DISTRIBUTION 

In a single-mode port cavity in the presence of losses, 
the scattering matrix S is a subunitary lxl matrix that 
we parametrize as 



S = VRe w . 



(2.1) 



Here R is the reflection coefficient and 9 is minus twice 
the phase shift due to the scattering process into the cav- 
ity. Those quantities are measured in actual experiments 
with microwaves. 

When the classical dynamics of the cavity is chaotic, 
S can be modeled by a random matrix distributed with 
a probability law. In the presence of both, absorption 
and imperfect coupling, the statistical probability distri- 
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bution of S is given by [If 



with 



(2.2) 



where /3 = 1 (2) denotes the presence (absence) of time 
reversal symmetry. Here, 



Kg® 



and (Rq) gives the probability distribution 
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1 - S(S) 
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of the lxl scattering matrix Sq = V^Rq el8 ° ■ This ma- 
trix describes the scattering of the system with losses but 
perfect coupling (see Fig. [I} . The parameters of So are 
the reflection coefficient Rq and the phase 9q , in an equiv- 
alent way as S. Note that 9o in Eq. (|2.4p is uniformly 
distributed between and 2-7T, while Ro is distributed 
according to p { o 0) (R o )- For (3 = 2 p]} 



x [ 7 ( e 7_i) + (i + 7 _ e 7)( 1 _ jRo) ] (2.5) 
while for any (3 (we extend the existing result for (3 = 1 



e -y/(l-Ro) 



Po (i?o) -(l— Ro? P ° {l^Ik 



fi(2){w) = / dt- 

Jw(0) 



■7*/2 



(1 +t) 3 / 2 



Si(2)H 
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(2.10) 



L(o) x/FH (i + 3/2 ' 
^1(2) H = / d< — /tM ■ 

j™(0) vt(l + i) 

x[7+(l-e-T)(7t-2)], 
/•oo(«i) 2 e~ 7 */ 2 

ji(2)H = / df , i TtTt 



Since Wi(a;) is very complicated, there are several at- 
tempts to interpolate p q 1 \Rq) between the two well 
known limits of strong (7—5-00) and weak (7 — > 0) ab- 
sorption 

The relation between S and So is given by the trans- 
formation fijl 



S(S ) 



(S) 



l + (S)S a ' 



(2.11) 



where (S) is the averge of S. In general (S), known as 
the optical scattering matrix, is a measure of the prompt 
responses in the system due to direct processes. In our 
case an imperfect coupling give rise to direct reflections. 
The coupling of the antenna to the cavity is quantified 
by 



T a = l-\(S)f 



(2.12) 



such that perfect coupling means (S) = or T a = 1, and 
5 reduces to So as well as pjg. (S) becomes the same as 

(2.6) Po 1 (So)- Then, by experimental measumerement of (S) 
we can calculate the coupling T a . 



where P (/3) (x), with x = (1 + Rq)/(1 - Ro), leads to the 
integrated probability distribution 

POO 

W {x) = / dxP^\x). 



(2.7) 



Wp(x) is the positive monotonically decaying function 



W 2 {x) = \e~ lxl2 [e 7/2 (x + 1) - e-~< /2 {x - 1)] (2. 



for (3 = 2 while for (3 = 1 
x + 1 



W 1 (x) 



4tt 



fi(w)g 2 (w) + f 2 {w)g 1 (w) 



+ h\(w)ji(w) + h.2(w)ji(w) 



w = (x-l)/2 



(2.9) 



III. THE ABSORPTION STRENGTH 7 AS A 
FUNCTION OF {R) AND T a 

The average (R) p can be calculated using directly Eq. 
()2.2p . However, it is necesary to write Ro as a func- 
tion of R and 9 (see Eq. (|2.3[) ), substitute the result- 
ing expression in the corresponding distribution p ^(Rq) 
above, and integrate with respect to R and 9. Because 
the resulting expression is still more complicated than 
,08) / 



p /S ,[(R), we prefer to calculate (R)/3 taking advantage of 



dP^(S) = dP^'(So) and integrate with respect to R Q 

and 60 ■ 

From Eq. (|2.4[) we have 



(R)p = J' T ^ j dR«R{R«Mp£\Ro) 



(3.1) 
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to experimentalists. We observe that to a single value of 
(R}/3 corresponds an infinite number of 7 values, each one 
corresponding to one value of T a . This means that the 
coupling has a similar effect as the losses on the scattering 
properties.^ e. the coupling mimic the absorption and 
viceversa [13j . This can be easily seen by writing the 
argument of the integral of Eq. (|3.3[) in powers of (1 — 
T a )Ro and performing the integral term by term. The 
result is a useful expression for (R)p, namely 



(R) fi = l-T a + T 2 a V(l - Tj n -V (1%) 



1 



where 



(3.5) 



FIG. 2: The average of R is a monotonically decreasing func- 
tion of the losses 7. We show 7 as a function of (R)p for 
several values of the coupling T a for (a) (3 — 2 and (b) (3 = 1. 
To a given value of (R)p it corresponds an infinite values of 
7, each for a T a given. 



where it remains to write R as a function of Rq and 9q. 
Using Eq. (T2~TTj) and that (S) = y/l - T a (the phase in 
the single mode case is not important) we arrive to 



R = 



R + (l- T a ) - 2 VT^TIVj^ cos g Q 
1 + (1 - T a )R Q - 2VT^7^ViuTcos O ' 



(3.2) 



If we subtitute this expression for R into Eq. (|3.ip and 
integrate with respect to 9q, the result is 



l-T n 



l-Ro 



1 - (1 - T a )R 



Po 0) (Ro)dR o , (3.3) 



where we use that plf\Ro) is normalized to unity. 

We can check the two limits of no coupling and per- 
fect coupling of the antenna to the cavity. For T a = 0, 
(R) ^ = 1 which is compatible with the argument that 
the waveguide is blocked and the wave never enters the 
cavity; hence there are no losses. On the opposite, T a = 1 
leads to {R)p = (Ro)p, which again is expected because 
the coupling is perfect such that S = Sq. 

When < T a < 1, the integral in Eq. (|3 . 3[) is not easy 
to do analytically but it can be done numerically. For 
(3 = 2 we directly substitute Eq. (|2"T5) into Eq. (573]) and 
integrate numerically with respect to Rq. The calcula- 
tion for /3 = 1 is more complicated. However, a simpler 
formula for (R)p for any (5 can be obtained by an inte- 
gration of Eq. (|3.3p by parts, giving the result (see App. 



{R)p = 1 - T a + 2Tl 



Wf3(x) 



1 [x + 1 - (1 - T a )(x - 1)Y 



rdx. 



(3.4) 

The remaining integrations can be done numerically. 

In Fig. [5] we show the results for (a) f3 = 2 and (b) 
(3=1 for several values of T a , as can be more useful 



Jo 



dR . 



(3.6) 



Again, T a — implies (R)p = 1, while (R)p = (Ro)p 
for T a = 1. The first two terms in Eq. (375]), 1 - T a , is 
the direct reflection due to the coupling at the entrance 
to the cavity, while the remaining terms are the reflec- 
tions after the multiple scattering occurs inside the cavity 
where losses are present. This becomes evident once we 
calculate below (R) p in both weak and strong absorption 
regimes. 

Then, in an actual experiment it is necessary to mea- 
sure both parameters, T a and 7. Once T a is calculated 
from Eq. (|2"32"j) . 7 is obtained from {R)p, Eqs. (pTB")) and 
(|3.4p , which is also obtained from the experimental data. 
Our formula allows to discern what part of the wave is 
reflected and what part is lost by absorption. 



A. Strong absorption regime 

In this limit, the probability density distribution of i?o 
reduces to H, fl7j | 



pW(fl ) = ^ r____ 



-apRo/0-Ro) 



(3.7) 



where we have introduced otp = 7/3/2. 

Substituting Eq. (|3.7|) into Eq. (13. 6|) it is easy to show 
that (see App. iBTjl 



(3.8) 



{R%)p » — — > as a/3 -> 00, 



which is consistent with the result (Ro)/3 = l/ a /9 °^ 
[T|. Then, Eq. ([33]) gives 



1-T„ 



(3.9) 



This means, the wave entering into the cavity with strong 
absorption never back and the only reflection occurs at 
the entrance to the cavity. 
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IV. CONCLUSIONS 

We have presented a semi- analytical formula to calcu- 
late the absorption strength 7, due to losses in a chaotic 
cavity, taking the coupling of the single channel port into 
account. Our result is useful in experiments with mi- 
crowave networks and one port cavities. We have shown 
that the imperfect coupling of the antenna to the chaotic 
systems has a similar effect as the absorption, and vicev- 
ersa, on the scattering properties. This formula responses 
to the necesity to calculate an accurate value of 7. We 
recall that the procedure one channel case can be applied 
for higher number of channels. 



FIG. 3: The weak and strong absorption limits of the average 
of R intersect at 7 = T a , fixing the criterion for low or high 

strength of absorption: (a) (3 = 2 and (b) (3 = 1. APPENDIX A: CALCULATION OF EQ. (f3^jl 



B. Weak absorption regime 



In this regime [171 ]. 



Po K H o) - 



a a +l3/2 e - ap /(i-ao) 



r(l+/3/2) (l-i? ) 2+/3 / 2 ' 



(3.10) 



where T(x) is the Gamma function [201 ]. Substituting into 
Eq. (|3.6p . after some simplifications, we get (see App. 
ED 



(3.11) 



2n 

that we insert into Eq. (|3.5[) to have 
T 2 T 



(R) 



1-T„- 



l-T c 

OO 

7T a 2 ^n(l-T a )- 1 



— y^(i-T a )» 

-L In „ 
n— 

(3.12) 



The fourth term on the right hand side is just the infi- 
nite geometric series and the fifth one its derivative with 
respect to (1 — T a ). They can be summed to give 



(R) « 1 - 7. 



(3.13) 



In this limit a small part of the wave that enters into the 
cavity is lost by absorption, and the reflection is lightly 
less than unity. 

Eqs. (|3.9|) and (|3.13|) represent the average of R in 
both limits of strong and weak absorption, respectively, 
where they arc independent of the symmetry (3. The 
intersection of both limits occurs at 7 = T a (see Fig. 
[3]). Then, the criterion for strong or weak absorption 
becomes 7 > T a or 7 < T a when the coupling is not 
perfect. 



We start with the substitution ofp^ (R Q ) of Eq. (|2.6|) 
into Eq. (|3.3|) . Making the appropriate change of vari- 
ables, the result is writen as 



where 



(R)p = l-2T a I , 



I f , = I A{x)P^\x)dx 



(Al) 



(A2) 



with Pq(x) giving rise to Wp(x) as in Eq. (|2.7p and 



A(x) 



1 



(A3) 



a; + l-(l-T )(a;-l) 

We integrate by parts identifying 

. / n dA(x) , , , 

u = A(x), du = — j^-dx, (A4) 
dx 

v = -Wp(x), dv = P^\x)dx. (A5) 
The result of the integration is 

I f3 = - A(x)Wf,(x)\™ =1 + ^ W (x)^-dx, (A6) 



For (3 = 2, Eqs. (JMJ) and (55) gives 



(A7) 



We show below that the same result is valid for (3=1. 

First, we evaluate A(x)W\{x) at x = 00, or w — 00, 
using Eqs. (|2.9p and (|A3|) . For instance, we consider 
the term fi(w)g2(w) appearing in Eq. (12. 9p . Defining 
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y = t/w we can write this term as 
f 1 (w)g 2 (w)\ 



(1-e-n lim r \ V{V \ } <""- lwy/2 
w ^°°Ji y (1 + wy) 

1 we ^wy/2 dy 



o y/y(l - y)(l + wy) 3 
lim / ,/— —we-W^y 



i yy(i + wy) 3 



(A8) 



Here, we make use of a definition of the Dirac delta func- 
tion, namely 



lim we^ wv/2 = ~5(y). 



(A9) 



The interval of integration in Eq. (|A8|) does not 
include the argument of S(y); as a consequence 
fi(w)g2(w)\ w=co = 0. In a similar way it can be shown 
that the remaining terms in Eq. (|2.9[) gives zero when 
evaluated at w = oo. Then, 



A(x)Wi(ir)| x=00 =0. 



(A10) 



Now, we consider the evaluation of A(x)Wi (x) at x — 
1, or w = 0. From Eqs. (|2. 10|1 it is easy to see that the 
first term in W\{x) gives 

h{w)g 2 {w)\ w=0 = lim / 7=== = „ cfy 



oJo v/y(l-y)(l + wy) 3 

dt, (All) 



which reduces to 

/•oo te -7*/2 / _ ]\ 

/iHa!HI«=o = ^ y o ( 1 + f ) 3 /2 ( 1 " e 7 + 7 J dL 

(A12) 

Similarly, we may also shown that the third term in 
Wi(x) gives 



hl(w)j 2 (w)\ w=0 = 7T 



-7*/2 



VT+t 
x [7 + (1 - e~ 7 ) ( 7 i - 2)] dt. (A13) 

The second and the fourth terms in W\{x) are zero 
when evaluated at x = 1. We consider the second one 
only. It is 



/ 2 (w)ffi (w)\ w=0 



lim 



\ ? {1 - V lv?e-™"> (l 



oJ y (1 + wy) 3 

/•oo e -7*/2 

L y/t(t - W)(l + t) 3 



1 



rft, 



wy) 
(A14) 



which reduces to zero. From Eqs. (|A12|) . (|A13|) . and (|A"3 

we get 



^)Wi(a:)U= 1 = i(l- 



-7*/2 



(l + 3/2 



- 7 t/2 



+ (7 - 2 + 2e~ 7 ) 



- 7 t/2 



/o (1 + t) 3 / 2 W J Jo (1 + t) 1 / 2 

/•<» fp-lt/2 

+ -' {1 - e ' Y} l (TTW (A15) 

Integrating by parts, we can establish the following 
relations 



■dt, (A16) 



00 e -7*/2 2 1 Z" 00 e~ 7 */ 2 

(1 + 1/2 ^ = 7~7/o (l + t) 3/ 

00 ie~ 7 */ 2 , /2\ 2 If 00 te~ 7 */ 2 , 

= / ~, ~TT77T dt 



o (l + t)Va-" ■ V7; 7 7o (1 + i) 3 / 2 
2 /- 00 e - 7 */ 2 



7 2 7o (1+i) 3 / 2 



dt. (A17) 



Substituting Eqs. (|A16j) and (|A17|) into Eq. (|A15j) we 

obtain 



A(x)Wi(x)\ 3 



(A18) 



Therefore, using Eqs. {£7|>, (|AT0)) . and (fATg)) . Eq. (|A6]l 

can be written as 

J^ + jTw^^dz, ( Al9 ) 

Finally, Eqs. ((AT]) . (|A3|) . and (|A19|) gives the result 
shown in Eq. 1|3.4I) . 



APPENDIX B: CALCULATION OF {S^) p 
1. Strong absorption limit 

In this limit, Eq. (13.71) can still be simplified to the 
Rayleigh distribution [3j, [l9| 



(Bl) 



which substituted into Eq. (|3.6[) lead us to 



a [ RSe- a ^ Ro dR 
Jo 

-e-^+—{B^-%, 

a/3 



(B2) 



where an integration by parts was done. This expression 
can be iterated to obtain the precise result 



n-l 



— ' (n — m)!a™ ' 



(B3) 



where the limit ap — > 00 gives the result shown in 
Ea. dm 
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2. Weak absorption limit 

We substitute Eq. (|3TT0]) into Eq. (J3T6J) to have 



1+0/2 

j3 



1 e -ap/{l-Ra) 
Rri ~ ^ s r, , a /o d Rq 



1+0/2 

<8 



r(i + /3/2) 



e -afl*(a; _ Ijnjff -"rf^ (B4) 



where we have change to the variable x = (1 — Rq)^ 1 . Us- 
ing the binomial expansion for (x — 1)™ we write Eq. (1B4|) 
as 

(B5) 



where T(a,x) is the incomplete Gamma function [20| . 
Keeping linear terms in only, we arrive to the result 
shown in Eq. (|3.1ip . We recall that special care for (3 = 2 
should be taken because of the divergence of the Gamma 
functions, however result in Eq. (|3.11[) is valid. 
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